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■ Abstract. We present new results and give a concise review of recent previous 

results on the asymptotics for large spin of the low- lying spectrum of the ferromag- 
netic XXZ Heisenberg chain with kink boundary conditions. Our main interest is 
to gain detailed information on the interface ground states of this model and the 
low-lying excitations above them. The new and most detailed results are obtained 
f*"^ I using a rigorous version of bosonization, which can be interpreted as a quantum 

' central limit theorem. 

o 
m 
o 

O^- 1. Introduction 

, In recent years the XXZ model has become a popular model to study properties 

of interface states in quantum lattice models. As an interpolation between the Ising 
model and the isotropic (XXX) Heisenberg ferromagnet, the ferromagnetic XXZ 
model has the interesting features of both. By considering the Ising model and the 
XXX model as limiting cases of the XXZ model, intuition about these two limits 
\ can be used to better understand the XXZ model. In this paper we are interested in 

the large-spin asymptotics of the low-lying excitation spectrum of the XXZ chain, in 
particular the excitations above the kink (or interface) ground states of the model. 
In a nutshell, our main result is that the spin-wave approximation, in the sense 
of Dyson |6j, becomes exact in the limit of infinitely large spin. The technical 
statements are given in Section |21 First, we introduce the model and the main 
notations and give a quick summary of the relevant previous results. 
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For J = 1/2, 1, 3/2,..., the spin- J XXZ Hamiltonian on an interval A = [a, h] C 
with kink boundary conditions, is given by 



b-l 



Hj^K = (1-1) 



x=a 



with 5* the spin- J matrices acting on site x: 

We will also use the spin raising and lowering operators: 5+ and , by = 

^x ^ ^^x- 

We begin with a brief overview of the main results obtained for the Hamiltonians 
Hj\. The spin 1/2 model, J = 1/2, is Bethe Ansatz solvable and posseses a quantum 
group symmetry ^H]- Consequently, there are a number of results specific to the 
spin 1/2 case (E.g., see Pl llU lTH]). Since the main focus in this paper is on large- J 
behavior, we will not discuss these specific results here. 

The Hamiltonian p.lj) is symmetric under global rotations around the 3-axis 
generated by = J^xeA'^xy which represents the third component of the total 

magnetization. Hence, Hjx is block diagonal, and it is known that in each sector 
corresponding to a given eigenvalue of S^^ot a there is exactly one ground state, i.e., 
in each sector is a simple eigenvalue [I^ : 

Hj,A<''' = 5L.A$f ) = M$f ) 

where M = — |A|J, — |A|J+ 1,...,|A|J. The (unnormalized) eigenvector is 
given by 

<'= E n C !■!„)' V-'K".}) 

{m,} xGA ^ 

where we introduced the parameter g, < g < 1, by the equation 2 A = q + q~^. 
These ground states have a magnetization profile that shows an interface, or kink, 
with a location depending on the value of M. For a short review on the properties 
of these ground states see . 

In many instances it is important to consider the thermodynamic limit, i.e, the 
limit of infinitely long chains. To this end, we consider a strictly increasing sequence 
of numbers a„ G Z'^ and lim„ a„ = oo, and a sequence of volumes A„ = [— a„ + 1, a„]. 
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The set of eigenvalues of 5'^^ot.A„ is then 

Mn = {-2a„ J, -2a„ J + 1, . . . , 2a„ J} 

and since 2 J is an integer, we have C A^m for n < m. Hence we can fix M G Z, 
take no large enough such that M E M.n for all n > tiq, and consider a sequence of 
states 



\^A„ ' ^A„ / 



It is shown in [TIH IT^ that in the limit n — > oo, any such sequence converges in 
norm to a unique state u^^^ on the quasi-local algebra of observables 21 which is the 
norm completion of the algebra of local observables given by 



2lioc = U ®xeAM2j+i(C) 



and each lu^^"^ is a ground state for the derivation 5j defined by 

6j{X)=limt[Hi,X] 



I.e., 

JM) > for X G 2lioc 

All these infinite volume kink states have the same GNS Hilbert space, namely 
the incomplete tensor product Hilbert space 



where 



f2y 



- J) if y < 
|J) ify>0 

Also denote fl = ®j;gz ^ 'Hj, and the GNS Hamiltonian on Hj will be denoted 
Hj. 

It was proved in p!3tl2T| that, for all J = 1/2,1,3/2,..., these Hamiltonians 
have a gap above the ground state eigenvalue, which is 0. Let us denote the gap by 
7j Af In the case of J = 1/2, the exact value of the gap was previously known to be 
1 — A~^, for all A > 1. In ^2] a very explicit conjecture is made about the value 
of the gap in the limit J — > cx). For all finite J, it is a periodic function of M, with 
period 2J. The conjecture in is as follows: 



Conjecture 1.1. For all fx e R, the limit 

lim ^-fj,^j = 7(r) 



4 



TOM MICHOEL AND BRUNO NACHTERGAELE 



exists and given by the smallest positive eigenvalue of the Jacobi operator h^'^\ defined 
below in (|2.1|) . where are is the solution of the equation 

+00 

/i = tanh(?7(a; — r)) 

x=—oo 

with 7] determined by A = cosh?]. 

A partial result towards this conjecture was proved in @|. Namely, there it is 
shown that there are constants ci > 0, and C2 > 0, independent of M and J, such 
that 

ClJ < 1J,M < C2J . 

In this paper we prove that the value of the gap claimed in the conjecture is 
asymptotically correct. The results in this paper by themselves, however, do not 
amount to a proof of the conjecture as stated. Roughly speaking, we obtain the 
result in the "grand-canonical ensemble" , and with the aid of a ground state selec- 
tion mechanism that localizes the kink. The conjecture is stated in the "canonical 
ensemble", i.e., for fixed fi = M/J, in which the kink is automatically localized at 
a fixed location. As is often the case the distinction between canonical and grand- 
canonical results seems merely technical, but proving mathematical equivalence of 
both formulations is often highly non-trivial. In fact, equivalence of ensembles in 
the usual sense does not hold in the present situation. To prove the conjecture as 
stated above, some additional work has to be done. We will report on this further 
work in a future publication 16 . 

The conjecture of fW\ was based on results form perturbation theory and nu- 
merical calculations on small systems presented in j2I], as well as on a heuristic 
calculation leading to a Boson model. 

The idea is to apply a rigorous version of Dyson's spin wave formalism to the XXZ 
chain. Mathematically speaking, the task is to control the quadratic approximation, 
described by a quasi-free system of Bosons, and show that this approximation be- 
comes exact in the limit J — 00. 

Several authors have attempted to do this for the XXX model, with interesting 
results pil22j. In these works, the authors considered the XXX model in an external 
magnetic field, and it was necessary to let the strength of the field diverge as J — > 
00. Such a field selects a particular ground state (out of the infinite number of 
them), and creates a gap in the spectrum. This allows one to proceed, but it 
limits the mathematical applicability of the spin-wave formalism. In the case of the 
XXZ-model, the situation is somewhat better. First, the XXZ chain by itself (i.e., 
without external magnetic field) already has a non-vanishing spectral gap. Second, 
although the infinite XXZ chain also has an infinite number of ground states — with 
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the degeneracy now corresponding to the arbitrary position of the kink — , any field 
at just one site with a non-vanishing component in the XY-plane, a so-called pinning 
field, will select a unique ground state, for finite J jHI- Moreover, the magnitude of 
this field, as we will show, can be taken of smaller order in J. These features of the 
XXZ model will allow us to prove asymptotic properties of the model itself. 



2. Main results 

2.1. The limiting Boson model. Our main result will be that the spectrum 
of the XXZ chain, in the limit of infinite spin, can be understood as the spectrum 
of a non-interacting system of Bosons on the chain, with one-particle Hamiltonian, 
U'-\ defined on fij) by 

(/iWt,)^ = e^v^ - ^(^x-i + v^+i) (2.1) 

where 

2 cosh(?7(x — r)Y 



cosh(?7(x — 1 — r)) cosh(r7(x + 1 — r)) 

with Tj = — Ing or, equivalently, A = cosh?], and r G IR is the position of the 
kink in the reference ground state, /i^'"^ has the form of the discrete Laplacian 
(kinetic energy) plus a diagonal term given by e^, which is an exponentially localized 
potential well centered around the interface. 

We list some properties of h^'^\ some of which are easily proved, while other 
more detailed properties about its spectrum have at this point only been verified 
numerically. We will discuss these in more detail elsewhere. 

(i) h'^'^hs a positive operator and hence the Fock state is a ground state for 

(ii) h^^^has an eigenvalue with eigenvector vq G /^(Z), up to normalization 
defined by 

Vqx = sin 9^J^ 



cosh.(ri{x — r)) 



(iii) the bottom of the continuous spectrum of h^^^ is given by 2(1 — A ^); 

(iv) the first excited state of /i'^'"^ corresponds to an isolated eigenvalue 7*^''^ below 
the continuum; 

(v) the (A"^, r)-plane is divided in a region where 7Wis the only eigenvalue be- 
low the continuum, and a region where there is another isolated eigenvalue 
between 7^'')and the bottom of the continuum. 
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Let h'^'^^ (x, y) be the bi-infinite Jacobi matrix expressing ()2.1|) in the standard 
Kronecker-delta basis of ^^(Z), i.e., 

The Boson Hamiltonian is then given by second quantization of h^'^^: 

x,y£Z 

where and ay are the creation and annihilation operators for a boson at site x 
and y, respectively. They act on the bosonic Fock space with one-particle space 
£^(Z), JF, and satisfy the canonical commutation relations (CCR) 

Let G denote the vacuum vector which, up to a scalar factor is uniquely 
characterized by the property a^^l = 0, for all a; G Z. 

We will often use the following standard orthonormal basis in JF. Introduce 

n = {n^ G N}^gz 

Then, the set {ipn \ n G Af}, where 

^^ = U7Au-M:r^ (2-2) 

is an orthonormal basis of JF. 

The GNS Hilbert spaces of the spin chains, TCj, J > 1/2,1,3/2,..., can be 
identified with a nested sequence of subspaces of JF, defined for each J, as the linear 
span of all vectors ipa, with <2J + 1. We will use this identification throughout 
the paper, and we will use the same symbol <^n to denote a vector in the spin Hilbert 
space Hj and the boson space JF. We will also use the projections Pn,x on JF which 
projects onto the first n boson states at site x, i.e., on the states ipn ()3.15|) with 
< ria; < n, and denote P„ = Hx-Pri,!:, i-e., Pj projects onto the states ipf^ with 
n G Aj, and hence 

Hj = PjJ^ 

More details on the Boson model are given in Section 13.31 
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2.2. Ground state selection by the external field. To prove full conver- 
gence of the low-energy spectrum, we need to add an external field. A physical 
external field would be of the form 

J2 ^J^- ■ (2-3) 

X 

For our purposes, however, the field is a perturbation and our results will generally 
be more interesting if we can proof them with smaller perturbations. Ideally, a 
vanishingly small fied localized at one site shoudl suffice to select a reference ground 
state. It turns out that we cannot quite do this in the present setup. We shall use 
a perturbation of the form 

s\xp{hj^x- S^) (2.4) 

X 

which is still significantly smaller than a uniform field. The meaning of this operator 
is clear if we express states in a basis which is diagonal for each hj^x ■ Sx- It is 
important to remark that we can take 

\hj^x\ = hj = h{J\YiJf'' 

such that after scaling with the external field vanishes, in contrast with the 
fields employed in previous works on the XXX model |21I2ZI- In fact some results 
can be obtained with hj = 0, or hj = hlnJ. 

Mathematically, the field ()2.4j) does slightly more than a field localized at one site 
can achieve. It not only pins the interface, but also puts some control on the local 
fiuctuations around the selected ground state, which cannot be controlled otherwise. 
The fact that we can let the field vanish as J increases, and that we do not need a 
global field like (j2.3p . are signs that these fiuctuations are smaller in the XXZ model 
than in the isotropic model. 

There is another mechanism of selecting a ground state, namely by restricting 
the full Hilbert space to a subspace of states with fixed total magnetization in the 
3-direction, since it is known that in each such sector there is exactly one ground 
state The ground states that are pinned by an external field are like grand- 
canonical averages of the canonical ground states with fixed magnetization. In the 
limit J oo the canonical description can be obtained from the grand canonical 
one by a result analogous to the result of [211 Section 5.11, 5.12] about equivalence 
of ensembles in the 2-dimensional, spin-^ XXZ model. Note however that there is 
no equivalence of ensembles in our situation, in the sense that correlation functions 
are typically different. But the difference between a canonical state and a grand- 
canonical state with the right average magnetization can be expressed completely in 
terms of the fluctuations of the total 3-magnetization in the grand canonical state 
which are non-zero even in the limit J — oo, while they are identically zero in 
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the canonical state. The results about the canonical description require significant 
additional work and will be discussed elsewhere [TB] . 

Finally we mention that because of the pinning field ()2.4|) . our results give only a 
partial proof of [121 Conjecture 2.5]. A full proof requires in addition that the lowest 
excited state can be obtained from the ground state by acting on it with 1 spin wave 
operator, or more generally by a finite number of operators independent of J. This 
is a problem that should be handled at the level of the spin system, rather than in 
the spin wave formalism. Some result in that direction have recently been obtained 
in [IHI. 

An advantage of the grand-canonical description is that it clearly exhibits how 
the boson limit arises as the first quantum correction to the classical limit. 

2.3. Statement of the main results. For A a self-adjoint operator and a,b E 
R, denote by P(^a,b){A) the spectral projection of A onto (a, b). For A acting on Fock 
space, denote 

• cr{A) the spectrum of A in JF; 

• CTj{A) the spectrum of PjAPj in Hj; 

• anj{A) the spectrum of PnjAPnj in Hnj, where Hnj is shorthand for Pn,,^- 

Also denote s-lim the strong, or strong resolvent, operator limit for bounded, resp. 
unbounded operators acting on 

In the GNS space Tij it is convenient to define S^^^ in the renormalized sense: 
Also denote 

= ^ [cos^^'^) - sgn(x - I)] 
i.e., /i is the 3-magnetization of a classical ground state {{Ox\fp)}xez (see section 

EH)- 

Theorem 2.1. We have 

s-lim hHj = H^'''^ 
srlim J Si, = fil 

J—>-oo 

The proof of this result is given in Section 

Corollary 2.2. (i) If X e a{H^'^^), there exists Xj e a,j{jHj) such that 

lim Xj = X 

J— >oo 
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(ii) Ifa,b E R, and a,b ^ app{H^^'^), then 

s-limP(,,,)(iifj) = P(,,,)(ijM) 

J — >oo 

(iii) Ifa,bE R, and /x G (a, 6), then 

S;limP(a,6)(j5,3ot) = 1 

This corollary is a direct consequence of Proposition 15.61 and \20\ Theorem 
VIII.24]. 

In addition we can obtain spectral concentration of jHj around discrete eigen- 
values of H^^^ (also proved in Section ing . 

Theorem 2.3. For every isolated eigenvalue E of H^''\ there exists an interval 

Ij = (E — £j, E + €j) with lim ej-^ = 

J^oo In J 

such that for any interval I around E s.t. ID a(H^'''^) = {E}: 

s-lim Pi\i ^ H j) =0 

Applying the same reasoning to S'^.t, we find that the interval (a, b) in item (iii) 
of Corollary 15.71 can be chosen as 

(a, 6) = (/i - ej, /i + ej) 

with again limj ejj(ln J)^^ = 0. 

To prove full convergence of the spectrum, we have to add the external field (|2.4|) 
to Hj, or, to have a positive operator, add: 

Fj = sup(|/ij,^.| J - hj^^ ■ S^) = hjsnp{J - S^) = hjsup{N^) (2.5) 

X XX 

with hj > 0, and A^^, = 0*0^.. Let us assume we add this field to jH-^ , so hj already 
contains the factor J~^. 

Take < nj < J as before, we get that on Tinj"*" H ?ij 

jHj + hj sup > hj sup > hjUjl 

X X 

Clearly, by choosing hj such that 

lim hjuj = 00 

J^oo 
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statements about the spectrum on Tij reduce to statements about the spectrum on 
Tinj- Or, if one chooses to make statements about the spectrum below a certain 
value E, it is sufficient to choose hj such that limj hjUj > E + e. 

Theorem 2.4. Let rij = [{J{\n jy'^f^']. If X ^ (t{H^'^), then X ^ anjijHj) for 
J large enough. 

This result is proved in Section 1331 

Hence we get convergence of the spectrum of jHj + hj sup^ A''^, if 

J \i/3 



lim hj = lim hj - — - = oo 

3. Derivation of the Boson limit 

3.1. The classical limit. The Boson limit can be considered as the ffist quan- 
tum correction the the classical limit. Therefore, we ffist discuss the classical limit. 

It is well known that for any quantum spin system, after rescaling each spin 
matrix by J^^ and taking the large spin limit J ^ oo, one obtains the corresponding 
classical spin system P^. For the XXZ chain in a finite volume this is defined by 
the Hamiltonian 

Ht{{a.U^)=Y,^--iay,^, + ay,^,)-alal^, + VT^A=^{al^ (3.1) 

x=a 

where is a unit vector in IR'^. Minimizing this function with respect to yields 
zero-energy configurations that are plane waves j2I], i.e., in spherical coordinates 
we find configurations ax\f) = {6x\<f), f G [— vr, vr] (the same at all sites), and 

eP = 2 arctan(g^-") 

where r G IR determines the value of the total 3-magnetization. 

Defining 77 = — Ing or A = cosh(?7), we have 

= TTfS4 = ^^Mvi^ - r)) (3.2) 

3i^g(^r) ^ ^g^-'' _ 1 ^^^^^ 

I _^ g2{x-r) cosh [rj {x — r)) 

such that the zero-energy solutions clearly describe kinks centered around r. 

For the classical model, to look at the low-energy behavior amounts to making 
a quadratic Taylor approximation to p.ip . At each site, the angle coordinates are 
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replaced by new coordinates 

qx = dx- Px = sin e^j'^ {ip^ - V?) 

and the resulting harmonic oscillator Hamiltonian is 

1 1 

= 2 + ~ -^{(Ixqx+i + PxPx+l] 



x=a 

b 



+ ^ Yl (^^ + ^^■) ~ ^(^^^^-1 + PxPx-l)) (3.4) 

x=a+l 

where are given by 



e 



j_ sin 6*^21 cosh(?7(x — r)) 



Asin^i*"^ Acosh(r/(x±l -r)) 
This can be derived using the identities of Lemma f5. II below. 

At sites other than the boundary sites we have a single-site potential 

_^ ^ _ 2 cosh(77(x — r))^ 

^ ^ ^ cosh(?7(x — 1 — r)) cosh(?7(x + 1 — r)) 

which is an exponentially localized well centered around the interface. 

3.2. Grand canonical states. For {6x, (fx) a general unit vector on the sphere 
at site X, we can define the coherent spin state in (C^'^"'"^)^ (see [HIT^): 

|(0.,^.)) = e^^^(^-'^^-^^"^"^^)|J) 
/ 2T \^^^ 

= E 7 (cosie,)-^+™^(sin|^,)-^-"^^e^('^-™^)'^^ |m,) 

This is particularly interesting if we choose the unit vectors at each site to be 
the classical zero-energy configurations. In a finite volume A, it is easy to see that 

® e-^^-|(e\¥')) = -^4^^ 

II^A II 

where \&^^is the generating vector for the ground state vectors ^\^\ i.e., the grand 
canonical ground state: 

M=-\A\J 

evaluated at z = q^'e'^'P = e-^'^e"^'^. 
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The fact that a classical ground state yields an exact quantum ground state 
through the coherent state representation, is because Hjj^ is a normal Hamiltonian 
in the sense of jH] , and the classical and quantum ground state energies are (exactly) 
related by the scaling factor .P. 

Since these states are product states, their thermodynamic limit is easily ob- 
tained. In the GNS Hilbert space Tij, define the embedding of |cr^^^ as 

xGZ\A 

Lemma 3.1. For a sequence of intervals A„ = [— a„ + l,an] tending to Z, we 
have for m > n 

, , , - 1 _ ^2(a,„-a„) 



Proof. 



xeAm\^n a;GAm\A„ 



=2- n fi 

= 2 - 2 JJ o2(x+r)y n (1 + o2(x-r)y 

<2(1 — e ^3^=^11 y z.a:=a„+i y 



x=a„+l 

1 — ^ 

where we used the inequalities (for "U > 0) and 1 — e~" < m. □ 



It follows that the sequence \o'^^Jk}p)) has a limit \o^~^\ip)) in that we can 
formally write as 
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In general, if we write z = \z\e^''''^, then 

SO it will be sufficient to restrict our detailed analysis to the grand canonical states 
^ ^ . The expectation in this state will be denoted uj^j[^ : 



^ II T fe-'J'') l|2 



and its thermodynamic limit 

In this case, the coherent states are rotations of the 'top' state | J) through an 

(r) 

angle 6x around the 2- axis. 

Introduce the notation = {9x \ 0), or in Cartesian coordinates 

«i^) = (sine\0,cose^) 

In the remainder, we will always keep r fixed and do not make explicit the 
dependence on r of various quantities. Notice that by periodicity it is sufficient to 
take r e [0,1). 

Denote by {e^, e^, e^} the standard basis in (the same at every site), and 

f = 

fl^ut^= sin 9^:^ e^ + cose^ el 



fl^ cos 9^:^ e^-sin^W el 



Conversely 



el^cos9^:^ fl + sm9^:^ fl 

^x fx 

el = -sm9^:^ f,+cos9P f! 

i.e., {fx^fxifx} iorm an orthonormal frame for R"^ and {fx, fx} orthonormal 
frame for the tangent plane IR^ to the unit sphere at u^x^. 

For Vx G R^, we denote by Vx G IR^ the projection of Vx onto the tangent plane 
at Ux^ (shifted to the origin), i.e., Vx = {vl,vl) and 



'i^l^Vx- fx = cos e^J^vl - sin O^vl ^ Vx ■ el ^ 



2 
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sl 




'si 


= fl- 


' Sx 


= cose^;^sl 




si 




'si 


= Sx 








si 




'si 


J X 


' Sx 


= sin si 


+ cos e^;^ si 



Conversely, if Vx G we associate to it a vector Vx G by putting the component 
along the wl^'^-axis zero: 

vl = cosO^^vl vl = vl vl = — sin 6'^'''' 5^ 

Also denote 

Vx ■ Sx = vlSl + vlSl + vlSl 
and define rotated spin operators (1, Eq. (3.9)]) 

(3.6) 
(3.7) 
(3.8) 

Hence we find that |(6'i^'',0)) is the 'top' state for the rotated spin operators: 

Sim\0)) = Sle-^''^'''^J) = e'^''^''''^Sl\J) = JKC^O)) (3.9) 
The rotated spin raising and lowering operators are 

Si: = Si ±iSl = - sin O^J^Sl + cos O^^^Sl ± iSl 

or 

St = - Sin e^:^Sl + cos' \- St - sin^ \- 5" 

n{r) n{r) 

S- = - sin el!-^Sl - sin^ ^ 5+ + cos^ ^ 5" 

One of the main observations is that because of ()3.9|) . it is much more convenient 
to introduce the spin wave formalism in the rotated spin basis than in the original 
one. Following ^2], introduce 



n = {ux e N}^.ez 

A/" = |n I n^. < oo| A(/ = |n I Vx: ra^ < 2J, ra^ < oo| (3.10) 



1 /2J^~5 



x£Z ^ \ 



The set {(ffi \ n E Afj} is an ortho normal basis for Tij. 

We conclude with a little lemma that complements ()3.5j) . 
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Lemma 3.2. For -n < < n, 

e'A,A\a^^)) = (^cos(i(^) + zcos^Wsin(i(^)) e"' ^^^^a^.M sinews.- j^W^ 

where 



Proof. We write the disentanglement relation Eq. (A4)] 

^iipSl ^ gi¥p{-isin0W{5++5-)+cosei'-)s|) ^ ^-iy-S- ^{lnyi)Sl ^iy+S+ 

where 

1/3 = (cos^v?) + ^cos6'i'■) sin^v?)) 

singj"^ sin(|y) . 

y+ = y- = — ; . = ax(<^) sm^^ > 

cos{^(f) + « COS sm{^(f) 

and recall that \(J^l^) is the product state of 'top' states for the S'-operators. □ 
It is now also clear how to choose the external field hj^^ in (j2.3j) such that 

namely 

hj,x = -hM^\ hj>0 . 

3.3. The boson chain. We consider immediately the infinite volume situation. 
Consider the Hilbert space of wave functions /^(Z) which we alternatively consider 
as the usual complex Hilbert space with inner product 



{v,w) = y^^vxt 



cWx 

xez 

or as a real linear space with symplectic form a and complex structure JT", i.e., 

venz) = {{vl,vl)eR')^^^ 

I \ 12 2 1 

The CCR -algebra CCR(/^(Z), a) is generated by unitaries {W{v) \ v ef{T)} which 
satisfy the commutation relations 

W{v)W{w) = e-^"(^'"')vr(t; + w) 
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The Fock state a) is the quasi-free state on CCR(/^(Z), a) determined by 

u{Wiv)) =e-5<''''') 

Its GNS representation is the usual Fock representation on a Fock space JF with a 
vacuum vector Q = ®xgz \0)x and creation and annihilation operators a| such that 

W{v) = e*^=^GzK<+^^"^) 

(We don't bother to distinguish between W{v) and its representative in the Fock 
representation). 

In this Fock representation we can define a quasi-free Boson Hamiltonian by 
canonically quantizing the classical harmonic oscillator Hamiltonian (|3.4p . This 
means replacing the position and momentum variables by canonical pairs Qx , Px, 
with commutation relations [qx,Py] = iSx,y and 

* _ Qx - iPx _ Qx + ipx 

The result is 

x=a x=a+l 

The corresponding infinite volume derivation is denoted ^'•^^(■) = limA/'z ^[-^a'\ ■], 
and the GNS Hamiltonian is denoted H^'^\ 



jj(r) ^ "^Sxalax - A ^a*(aa;-i + ax+i) 



X&i 



Denote a*{v) = '^x'^^^x- If is local, i.e., has only finitely many Vx non-zero, 
then 



where W^vs the bi-infinite Jacobi matrix defined on /^(Z) by 

{U''\)x = SxVx - + Vx+l) (3.12) 

For the finite system localized in A = [a, b], we have i[H^\ '^*('^)] = (^*(h\^v), where 
{h^Av)x = (h''^^v)x for X in the bulk [a + 1, 6 — 1], and at the boundary sites we get 

h\'Vb = £5 Vb - 
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Some important properties of h'^^^weie given in Section ITTl Recall the existence 
of a zero mode (property (ii)) given by: 

Vo,x = sin e^j'^ = — — ^ -- 

cosh[ri[x — r)) 

That this is an eigenvector of h^'^^with. eigenvalue zero, is easily verified: 

± 1 sin6'^2i • air) 1 • air) n 

^ A sm e^x' ^ 

The origin of the zero-mode vq is well understood. It arises from the rotation 
symmetry of Hj\. This can be seen from Lemma 13.21 by taking ip oc J~'^^^ and 
formally identifying J~^^^S~ with a* (this identification will be made more precise 
below). For every G N, there is a zero-energy vector tpQ^N corresponding to an 
iV-particle occupation of Vq: 

ipo^N = - — , J a*{vo)^n 



\vo\^VN\ 

We define Pq the projection on vq and Pq the projection onto the zero-energy vectors, 
i.e., 

Po = ^ (3.13) 

Po= ® |^o,iv)(^o,7v| (3.14) 

For completeness we mention that we will always use the standard orthonormal 
basis in JF, i.e., the set {(ffi \ n G A/"}, where now 



^ ^ [n 

x&Z ^ 



IIt^TTwiK^^ (3-15) 



We use the same symbol (pfi to denote a vector in the spin Hilbert space Tij and 
the boson space since we will use the identification of Tij with a subspace of JF 
as discussed before. 



4. The large spin limit as a quantum central limit 



Since the grand canonical states can be written as the 'all -|-' state for rotated 
spin operators, we are in the usual situation of a fully ferromagnetic state in which 
we expect a boson limit after rescaling with J"^''^, i.e., define 




in some sense. 
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One way to make this precise is to define fiuctuation operators: for Vx E R^, 
Fj{vx) = ^{v, ■ - u'^^'Xvx ■ S,)) 

i.e., Fjivx) measures the deviation from the ground state expectation value of the 
spin in the direction. Similar fiuctuation operators are used to study fiuctua- 
tions of extensive observables, and their thermodynamic limit can be taken as a 
noncommutative central limit 012] ■ 

A connection between the spin limit J ^ oo and these quantum central limits 
was made in pL^j, with the caveat that each spin- J had to be represented as a 
sum of spin- 1 's, instead of working with an irreducible representation. This latter 
restriction is however not necessary. We have, using results from P^ll4j: 

o;M(e--^^) = {cos(^|^;.|)+^^^^|^sin(l|^^ (4.1) 

uj^''\Fj{vx)Fj{wx)) = ]^{vx-Wx- ■ up){wx ■ uP) + i{vx X Wx) ■ uP) 
The latter quantity defines a (degenerate) inner product on R^: 

{vx, Wx)x = {Fjivx)Fj{wx)) (4.2) 
It is not hard to use ()4.H1 to show that 

lim J'Ve'^'^^'^A = e-s^'^-^-) (4.3) 

J^oo 

Clearly if either Vx or Wx is along the Mi^''-direction, then {vx,Wx)x = 0. Hence 
()4.2|) defines an inner product in [R^ = C, the tangent plane to the unit sphere at 
uP . If for Vx,Wx G [R^, Vx,Wx are the corresponding vectors in R^ (see Section 



{vx,Wx)x = {vx,Wx)x = {vl + ivDiwl + iwD 

i.e., the standard inner product in C. We see that there are no fiuctuations in the 
direction perpendicular to the tangent plane at the classical zero-energy solution. 
Two vectors in R^ at the same site will be called equivalent if their projection onto 

(r) 

the tangent plane at Ux is the same. 

For V = {Vx € IR^)xeZ) with only finitely many Vx 7^ 0, we simply extend this by 
putting 

Fjiv) = J2Fjivx) 

{v,w) = 2u[Fj{v)Fj{w)) = '^{vx,Wx)x 

xez 
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Using ()4.3j) . and standard techniques, the quantum central hmit theorem follows: 

lim a;('-)(e'^^("i) . . . e'^^("")) = u{W{vi) . . . W{vn)) 



J— »oo 



where a) is the Fock state on the CCR -algebra CCR(/^(Z), a) introduced in section 
13.31 and the vectors f i, . . . , on the r.h.s. mean their respective equivalence classes 
in /^(Z). The result can be understood intuitively from 



Vx ■ Sx — v^Sx + v^Sx — — iVx)Sx + + iVx)Sx 

for Vx G [R^ and corresponding Vx G R^. 

When studying properties of the GNS Hamiltonian, it is actually easier to make 
a correspondance between the GNS Hilbert spaces Hj of 10^"^^ and of uj. 

Following [22^, introduce the projection Pj^x on which projects onto the first 
2 J + 1 boson states at site x, i.e., on the states ()3.15j) with < Ux < 2 J, and 
denote Pj = Ylx^J,^^ projects onto the states ipfi with n G Aj, see (j3.1U|) . 

By identifying ipfi ()3.1H) with ipfi ()3.15p . it is clear that 

Hj = PjJ^ 

where = means unitarily equivalent. Under this equivalence, we find that the spin 
operators are given by P^ : 

^ Sx = Pja*xgj{x)^ ^='S'x = g.j{x)^'axPj J - Sl = PjaluxPj (4.4) 



^2J V2J 
where 



gjix) = gj{alax) 



and 



I - n<2J 
n>2J 



5. The low energy spectrum 

5.1. Some estimates for the Hamiltonian. We first need the following iden- 
tities: 
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Lemma 5.1. 

cos e^^l^ + cos ^iji = cos 9^^^ (5.1) 

A-i(sin^i'2i + sind^^l^) = e^sinO^;^ (5.2) 

sin e^;^ sin + cos 6^^ cos = (5-3) 

cos ^i") sin - sin ^^'^^ cos ^^i = tVI - A-2sin ^(j-) 

(5.4) 

^ sin sin + cos ^^i cos 9^;^ T Vl - A-2 cos ^^^^ = (5.5) 

Proof. These are straightforward computations using the definitions flH.2|l and 
(jH.Hp of cos e^J^ and sin6'i'^\ and the addition laws for sinh and cosh. □ 



With this lemma we can write the Hamiltonian in terms of the ^-operators. 
Corollary 5.2. 



x,x+l 



1 

2A' 



+ jVi^A^{cos9^:^SI - cos 5^+0 

+ v/r^A^(sine^5^^.\i - sme^:lslSl^,) 



where 



7.,.+i = et + v^r^A^cos^^ = e-^, - v^T^A^cos^it 



Proof. This follows immediately from the relations ()3.6p - ()3.8|) . and the pre- 
vious lemma. □ 

With the Hamiltonian in terms of the ^-operators, we can apply the unitary 
transformation (j4.4j) to write the spin Hamiltonian as an operator on JF: 

^ x=a 

b 



x=a+l 



+ 



v^r^Aj 

2J~" 



6-1 



cos^^^'-^iV^^-cos^WiVfl +5^7.,.+i 



2J 



+ 
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-2 



Pj (5.6) 



x=a 

and likewise for the oo-volume GNS Hamiltonian: 



-1 + gjix + iy'^a-^+i) 



lx,x+l- 



2J 



a;GZ 



P.J 



(5.7) 



In the usual language of spin wave theory [H | l7 | l22j . the first term in the Hamilton- 
ian is called the kinematical interaction -ffj^kin, and the second term the dynamical 
interaction Hj^dyn- The last term, which describes transitions between subspaces 
with constant number of particles, is not usually present. We denote it -ffj,tran. 
Note that we define these three operators with the right scaling already included, 
i.e., 

jHj = Pj {-ffj,kin + Hj^dyn + -f^J.tran} Pj 

We will only let these operators act on vectors in Hj, hence we may forget about 
the Pj. To further simplify some notation, introduce 

• the column vector A of annihilation operators: 



A 



ax 



• the diagonal matrix Gj: 

Gjix,y) = gj{x)5x^y 

Then we can write 

i^j,kin = A*GTW^GTA - ^ 

X 
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where h^"^^ is the one-particle boson Hamiltonian, see ()3.12|) . i.e., the matrix with 
entries 

In the following we will fix for every J an nj G N with < rij < J and make 
statements about the subspace Pnj'Hj of Hj. In the end we will formulate results 
on the whole of Tij by adding an external pinning field which will take care of the 
states in (PnjHj)^ fl Hj. For simplicity denote Hnj = Pnj'Hj. 

Lemma 5.3. On Tinj we have the lower hound 

where 7*-^-' is the spectral gap of h^^^ and Pq is defined in ()3.13j) . An upper bound 
(on the whole Tij) is given by 

Proof. projects onto the vectors with at most 2nj particles per site, such 
that on Hnj- 

Gj > gj{2nj)l 

Obviously Gj < 1 on Tij. The Lemma follows from the bounds on /i*^^^: 

and also Ex < 2. □ 

We are going to compare the spectrum of -ffj,kin with the spectrum of H^^\ Both 
operators commute with A^tot so it is sufficient to compare them on eigenstates of 
A^tot- Also for -ffj.dyn it is sufficient to look at eigenstates of iVtot- 

In the proof of the following lemmata we will use the following notation: for 
n, m G A'': 

m^ = nx + l 
T^n = m iff < m^^i = n^^i — 1 

my = Uy Wy X, X + 1 

^ _i ^ ,™ I rrix ~~ ^x 1 

n = m m ' 



X 



ruy = Uy Wy ^ X 
Lemma 5.4. Let i/jn e Hnj, \\^n\\ = 1; NtotipN = Nipj^. Then 
||(i^('^)-i^j,kin)^;v||<2(l + A-i)^ 
11^7 / II ^ 4njiV 

J.dynV^Af II S J 
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Proof. We can write 

= ^Cfi ipH with ^ \Cfif = 1 
n n 

where the sum runs over n e jVj for which n^; — N. On basis vectors, we have 

Hj,]^in(Pn = ^£xnxgj{nx)(Pn 

X 

- A"^ J^fK + l)9j{n:,)gj{nx-i - l)n^_i] '^Vr-n 

X 
X 

and for H^'^'> the same with the gj = 1. We compare term by term. The first one 
gives: 



X ) V'ri 



2 4iV2^2 



< 



J2 



where we used Sx < 2, Ux < 2nj and Ux — N . For the second term we use 
analogously gj{nx) > gj{2nj), and find: 

II J^c^^(|[(n^ + - [{rix + l)gj{nx)gj{nx-i - l)n^_i] '^^^</7j,-^ 



= ||X]^T^tim(["^^("^^-i + 1)]'''^ ~ [mxgj{mx - l)gj{mx-i){mx-i + l)]'''')v': 

m,x 

< |cT+_iml| + 1)]'^' - [mxgj{mx - l)gj{mx-i){mx-i + 

m X 



2;^E(Ei-T.t.JV.-. + i))(E 



and the same for the third term. Summing everything together we find 



^^'■^-i^j,kin)^7v|| <2(1 + A 



J 



For i?j,dyn we use the same reasoning, ^x,x+i < 1, and 

{ux - Ux+if < 2nl + 2n^+i < 4:nj{nx + Ux+i) 

to find 



injN 

J 
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□ 

For -ffj,tran we have the following estimate. 
Lemma 5.5. Lettp e Ttnj, = 1- Then 

where \vo\i is the l^-norm of the zero-mode. 
Proof. Let 

V' = XI ^^^^ 

n 

Then 

X 

n,x 

I 2 
n a; 

n X a; 

n X 

= \vo\i{2nj + l){Anjf 
and likewise for the second term. □ 

5.2. Strong convergence and spectral concentration. Recall the following 
definitions. For A a self-adjoint operator and a,b E R, denote by P(^a,b){A) the 
spectral projection of A onto (a, b). For A acting on Fock space, denote 

• (J {A) the spectrum of A in J^; 

• crj{A) the spectrum of PjAPj in T^j; 

• (7„j(A) the spectrum of PnjAP^j in Hnj- 

Also denote s-lim the strong, or strong resolvent, operator limit for bounded, resp. 
unbounded operators acting on J^. 
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In the GNS space Tij it is convenient to define Sl^^ in the renormahzed sense: 

a;GZ 

Also denote 

i.e., /i is the 3-magnetization of a classical ground state v^)}xgz (see section 

EH). 

Proposition 5.6. We have 

s-lim \Rj = H'^'^^ 
s-lim ^Sl^ = /il 

J—^oo 

Proof. Introduce the set T>, the finite linear space of vectors ipa with n G A/". 
D is a common core for jHj, for all J, and H^'^\ 

Take ip &T> arbitrary (but normalized for simplicity) and denote 

^ = 5Z ^^^^ 
ft 

= sup ^ n^. 

Note that by assumption, ip is a. finite sum of with Ylx''^^ ^ ^^"^ hence also 
< oo. 

Now take J large enough such that ip G From the proof of Lemma (5.41 and 
15.51 it is clear that we can use 2nj < N^, as soon as 2 J > N, hence 

11(1//. - H'M] < (3 + A-)f + 2VT^\,M,(i^f*±^p^y" 

The first result follows from fl^ Theorem VIII. 25 (a)]. 
To prove the second statement, write 

53 = cos 9^;^ Si -sin 9^;^ Si 

= Pj[cos9^:\j - AT,) - sin^M((7j(x)^a, + ^(^^(x)^) }Pj 
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X X 

Clearly cos6i^''Nx < Ntot, and for G as before 



sin ^(^'^)<c/j(x) 5 CH^fi) 



= ||^sin^(^'')c^(n^ + l)'/'^j(n^) 2(^a,^ 

n,x 

< {N^ + l)\vo\l 



Hence 



{jSL - <^{N^ + 2(iV^ + l)^^>o|i) 



and j5'|ot — ^ /^l strongly on JF. 
Corollary 5.7. 

(i) If \ E a{H^'''^), there exists Xj E aj{jHj) such that 

lim A 7 = A 



□ 



(ii) If a,b E R, and a,b ^ app{H^'^^), then 

J — >oo 

(iii) Ifa,bE R, and /x E {a,b), then 

Proof. The previous proposition and |2()| Theorem VIII. 24]. 



□ 



In addition we can prove spectral concentration of jHj around discrete eigen- 
values of H^''\ 

Proposition 5.8. For every isolated eigenvalue E of H^^\ there exists an in- 
terval 



Ij = {E — Ej, E + Ej) with lim £j 



J 



j^oo In J 

such that for any interval I around E s.t. ID a(H^'^^) = {E}: 

s-\imPj\j^{^Hj)=0 

s-limP,,(iff,)=P|^}(#«) 

J — :*00 
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Proof. Let iI)e be the simultaneous eigenvector of i^^*"^ with eigenvalue and 
of A'^tot with eigenvalue Ne- It follows that 

We have from Lemma f5 .41 

||(lH,_^;).fe|l<£(^ 

or, 

lim = 

J^oo In J " " 

Hence it follows that ii^ is a first order pseudo-eigenvalue with first-order pseudo 
eigenvector ipEi and the result follows from prP, Theorem XII. 22]. □ 

Remark 5.9. Applying the same reasoning to S^^^, we find that the interval 
(a, h) in item (iii) of Corollary 15.71 can be chosen as 

(a, 6) = {ji- ej,n + Ej) 

with again limj ejj(ln J)~^ = 0. 

5.3. Convergence of the spectrum with a pinning field. To prove full 
convergence of the spectrum, we have to add the external field fl2.4|) to Hj^ or, to 
have a positive operator, add: 

Fj = sup(|/ij,i.| J - hj^^ ■ Sec) = hjsup{J - S^) = hjsup{N^) (5.8) 



with hj > 0, see also the end of section I3T21 and = a^ttx- Let us assume we add 
this field to ■jH'^ , so hj already contains the factor J~^. 

Take < rij < J as before, we get that on Hnj'^ H Hj 

jHj + hj sup > hj sup > hjujl 

X X 

Clearly, by choosing hj such that 

lim hjUj = oo 

J— >oo 

statements about the spectrum on Hj reduce to statements about the spectrum on 
Hnj- Or, if one chooses to make statements about the spectrum below a certain 
value E, it is sufficient to choose hj such that limj hjUj > E + e. 

Convergence of the spectrum of ifj,kin + -f^j.dyn can be proved under the weakest 
assumptions on rij. First we prove convergence of the spectrum of Hj]^[^. 

Proposition 5.10. Let nj = [J(lnJ)~-^], where [■] denotes the integer part. If 
X ^ a{W^'^), then X ^ (Jnj{Hj\^i^ for J large enough. 
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Proof. Assume that A G crnj (-f^j.kin) for all J larger than some Jq. Take 5 > 
arbitrary, and ipj G Tinj a normalized approximate eigenvector: 

\\{Hj,^^-\)i,j\\<5 

More precisely, take iIjj such that 

-P(A-<5,A+5)(-f^J,kin)^J = V'J 

Since [i?j,kin, -^tot] = we can take ijjj an eigenstate of A^tot with eigenvalue Nj. 
Since tfjj is orthogonal to the ground state space, it follows from Lemma 15.31 that 

By our choice of rij, we have YmvjnjJ ^ = 0, or Ymvj gj[2nj) = 1, and for J large 
enough, 

-^<-^^ 

Putting this into the bounds of Lemma 15. 4t we find 

- X)^j\\ < - ffj,kin)V^j|| + ||(i^J,kin - 

and it follows that -ipj is an approximate eigenvector for H^^^ as well and A G 
ct(^W). □ 

Now we add Hj^yn- 

Proposition 5.11. Let again nj = [J(lnJ)"^]. If X ^ a{H^''^), then X ^ 
c^njl-f^j.kin + -f^J,dyn) foT J large enough. 

Proof. Since Hj ^yn > it follows that an approximate eigenvector ipj for 
Hj,],in + Hj,dyn Hiust Satisfy 

||-f^J,kin^j|| < X + 6 

with the same notation as in the previous proposition. Hence we get the same 
estimate on Nj as before, but this implies by Lemma (5.41 that limj || ifj^ dynV'jII = 
and tpj is an approximate eigenvector for H^^^ as well. □ 

Alternatively, these propositions prove that the spectrum of ifj^kin + Hj,dy^ + 
hjsnp^Nx converges to the spectrum of H^'^\ provided 

lim hj = lim = oo 

J J In J 

To add Hjtnm we clearly have to relax our condition on nj. 
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Proposition 5.12. Letnj = [{J{\nJ)-^y^']. If\ i thenX^ (^nA^Hj) 

for J large enough. 

Proof. By Lemma [5.51 we have for all if) G "Hnj 

and by assumption the r.h.s. goes to as J ^ oo. □ 

Hence we get convergence of the spectrum of jHj + hj sup^ if 

/ J \i/3 

lim hj = lim hj ) = oo 

The results obtained in this section include the statements made in Section IT^ 
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